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Abstract 

A general scheme is presented for controlling quantum systems using evolution driven by non- 
selective von Neumann measurements, with or without an additional tailored electromagnetic field. 
As an example, a 2-level quantum system controlled by non-selective quantum measurements is 
considered. The control goal is to find optimal system observables such that consecutive non- 
selective measurement of these observables transforms the system from a given initial state into 
a state which maximizes the expected value of a target operator (the objective). A complete 
analytical solution is found including explicit expressions for the optimal measured observables 
and for the maximal objective value given any target operator, any initial system density matrix, 
and any number of measurements. As an illustration, upper bounds on measurement-induced 
population transfer between the ground and the excited states for any number of measurements 
are found. The anti-Zeno effect is recovered in the limit of an infinite number of measurements. In 
this limit the system becomes completely controllable. The results establish the degree of control 
attainable by a finite number of measurements. 



Electronic address: apechen@princeton.edu 



'''Electronic address: hrabitz@princet on.edu 



I. INTRODUCTION 



A common goal in quantum control is to maximize the expected value of a given target 
operator through application of an external action to the system. Often such an exter- 
nal action is realized by a suitable tailored coherent control field, which steers the system 
from the initial state to a target state maximizing the expected value of the target oper- 
ator [l, [2I, li, 0]. A coherent field allows for controlled Hamiltonian evolution of the 
system. Another form of action on the system could be realized by tailoring the environment 
to induce control through non-unitary system dynamics p|. In this approach the suitably 
optimized, generally non-equilibrium and time dependent distribution function of the envi- 
ronment (e.g., incoherent radiation or a gas of electrons, atoms or molecules) is used as a 
control. Combining this type of incoherent control by the environment with a tailored co- 
herent control field allows for manipulation of both the Hamiltonian and dissipative aspects 
of the system dynamics. 

Quantum measurements can also be used as an external action to drive the system evolu- 
tion towards the desired control goal. There are two general types of quantum measurements: 
instantaneous von Neumann measurements (selective and non-selective) js] and continuous 
measurements 9J. If the measured operator is Q = ^2 i qiP-i, where qi is an eigenvalue of 
Q with the corresponding projector Pj, then the result qi of an instantaneous von Neu- 
mann measurement of Q is obtained with probability Pi = Tr [PiPs], where p$ is the state 
of the system just before the measurement. The state of the system just after the selective 
measurement with the result will be PjpsP/Pi- If a non-selective measurement of Q is 
performed (i.e., if the particular measurement result is not selected) then the system state 
just after the measurement will be ppsP- 

Non-unitary dynamics induced by measurement-driven quantum evolution was used re- 



cently in 



10| for mapping an unknown mixed quantum state onto a known target pure 



state. This goal was achieved with the help of sequential selective measurements of two 
non-commuting observables. After each measurement the outcome was observed and used 
to decide either to perform the next measurement or to stop the process. The same problem 
was studied in the presence of decoherence introduced by the environment ll|. The control 
of the population branching ratio between two degenerate states by continuous measure- 
ments was considered 12], while the effect of non-optimized measurements on control by 
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lasers was investigated [13]. 



Quantum measurements may be used in a feedback control scenario 14|, [l5|, [l6|, lljj]. In 
this approach continuous observations are performed, a controller processes the results of 
the measurements, and then based on these results modifications are made in the coherent 
control field in real time to alter the behavior of the quantum system. Optimal measurements 



19| , where the system state depends 



may also be used for quantum parameter estimation 18, 
on certain c- number parameters Q{. The goal is to find an optimal measurement strategy to 
extract the information about these parameters. 

In this paper we explore nonselective von Neumann measurements to control quantum 
dynamics. Any measurement performed on the system during its evolution has an influence 
on the dynamics. In particular, non-selective measurement of an observable with a non- 
degenerate spectrum acts on a quantum system by transforming its density matrix into 
diagonal form in the basis of eigenvectors of the operator corresponding to the observed 
quantity. Measuring different observables may produce different changes in the system's 
state. We optimize the measured observables such that their consecutive measurement 
modifies the density matrix to maximize the objective. The general formulation includes 
the use of optimal measurements along with a tailored coherent control field. A particular 
case corresponds to control only by measurements such that the coherent control field is 
not applied. For this case a complete analytical solution is found for a two-level system. 
Arbitrary target operators and initial states of the system are considered. The solution 
includes explicit expressions for the optimal measured observables and the maximal objective 
values attained. While control by measurements admits an explicit analytical solution for 
two-level systems, the generalization to the multilevel case is not straightforward. The 
situation becomes even more complicated if coherent control fields are used in addition 
to optimized measurements. For this case, numerical simulations are performed in 
for several quantum systems controlled by a tailored coherent control field together with 
optimization by learning control of quantum measurements. 

The quantum anti-Zeno effect 21[ can be used to steer the system from an initial to a tar- 
get state. In the anti-Zeno effect continuously measuring the projector E(t) = U(t)EW(t) 
steers the system into the state E(t), where E is a projector leaving the initial state un- 
changed and U(t) is a unitary operator. Continuous measurements in the quantum anti-Zeno 
effect are obtained as the limit of infinitely frequent von Neumann measurements. With the 
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anti-Zeno effect the system becomes completely controllable in this limit. In the laboratory 
it may be difficult to perform a large number of measurements in a short time interval. Thus, 
a balance may need to be struck between the number of measurements and the desired de- 
gree of control. In this paper we analytically establish the degree to which the system can be 
controlled by any given finite number of measurements. This result allows for determining 
the optimal control yield in balance with the cost of performing the measurements. 

The paper is organized as follows. In Sec. II the general concept of control by measure- 
ments is outlined. Section III presents the complete analytical solution for the problem of 
control by measurements in a two-level system, and as an application Sec. IV presents upper 
bounds on population transfer by non-selective measurements. In Sec. V the relation of this 
analysis with the quantum anti-Zeno effect is established. Brief conclusions are presented 
in Sec. VI. 

II. FORMULATION OF CONTROL BY MEASUREMENTS 

The control "parameters" in the present work are the observable operators Qi, • • ■ , Qn 
to be measured. The number of measurements N can also be optimized if the cost of each 
measurement is given. The scheme described here entails the consecutive laboratory mea- 
surement of the observables Q±, . . . , Qn on the same physical system, but the measurement 
results are not recorded and not used for feedback. The latter two restrictions could be 
lifted, if desired. 

Consider the effect of a non-selective measurement of an observable Q on the system 
density matrix. Let Q = qiPi be the spectral decomposition of the observable, where ^ 
is an eigenvalue and Pj is the corresponding projector such that P- = Pi, PiPj = 5ijPi and 

Pi = I- A non-selective measurement of the observable Q transforms the system density 
matrix p into A4q(p) := ^2 ii P i pP i . In particular, measuring an observable Q which has 
a non-degenerate spectrum diagonalizes the density matrix in the basis of Q. In this case 
■Mq(p) = J^PiPi, where p; = Tr [pPi] is the probability to get the value as the outcome 
of the measurement. 

There are classes of equivalent observables where measuring an observable Q makes the 
same transformation of the system density matrix as measuring any other observable Q from 
the same equivalence class. Two observables Q and Q are measurement- equivalent if for any 
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density matrix p one has M.q(p) = Aig(p). The observables Q and Q are measurement- 
equivalent if their spectral decompositions have the form Q = £\ gjPj and Q = g^p 
where for z 7^ j one has g» 7^ g.,- and g~j 7^ q~j. In particular, all observables of the form Q = gl, 
where I is the identity operator and q is a real number, are equivalent to Q — I. The latter 
observables are trivial in the sense that measuring any such observable does not change the 
system density matrix. 

Let po be the initial system density matrix. Consecutively measured observables 
Qi, . . . , Qn modify the initial system density matrix p into 

PN = M Qn o Mq n _^ ■ ■ ■ o M Ql (p ) (1) 

The typical goal in quantum control is to maximize the expectation value of a target operator 
assuming that initially the system is in a state pQ. The objective functional has the form 

MQi,...,QN]=Tr\p N e] (2) 

where p^ is defined by (JTJ). The control goal is to find, for given po and B, optimal observables 
Q° pt , . . . , Qtv* which maximize the objective functional to produce J™ ax . 

The general case also includes a tailored coherent electromagnetic field e(t) as a control 
where the dynamics of the system is governed by the two forms of external action: (a) mea- 
surements of observables Qi, ■ ■ ■ ,Qn at the times t%, . . . , tjsr, respectively, and (b) coherent 
evolution with a control field between the measurements. The former action induces non- 
unitary dynamics in the system. The latter action produces unitary evolution of the density 
matrix p(t) between the measurements according to the equation 

^ = -i[tf -^(t),p(t)] (3) 

where H Q is the free Hamiltonian of the system and p its dipole moment. The solution 
of Eq. Q at a time t, with an initial condition p(£o) = p &t t < t is given by a unitary 
transformation of the initial density matrix denoted as t/[t ,t)(p)- I n this notation the system 
density matrix at a target time T > t jv after N measurements will be 

P( T ) = u \t N ,T) o M Qn o U [tN _ lttN ) o M Qn _, o • • • o U [tut2 ) o M Ql o U [0}tl) (p ) (4) 

The density matrix p(T), which is dependent on the control e(t) and Qi, ■ ■ - Qn, determines 
the objective functional of the form J[s(t), Qi, . . . ,Qn] = Tr [p(T)9] with some target ob- 
servable 9. Here, in addition to the coherent field e(t), the observables Qi, ■ ■ ■ ,Qn are 



included as variables to be optimized. This general case is difficult to treat analytically, and 
for some models numerical simulations may be performed [^(J. In the next section we show 
that control only by measurements admits an analytical solution in the case of a two-level 
system. 

For an atomic multi-level system, practical measurements of the energy level populations 
can be performed using coherent radiation. For example, measuring the population of energy 
levels |1) and |2) of a two-level system can be performed by coupling the ground or the excited 
level by a laser pulse to some ancilla upper level and then observing the spontaneous emission 
from the ancilla level to a lower energy level. Such a measurement is described by projectors 
Pi = |1)(1| and P2 = |2)(2| and corresponds to measuring an observable of the form Q = 
q\P\ + g2-p2 with q x 7^ q 2 . This case shows the distinction between the use of coherent 
radiation for control and for measurements. The coherent radiation used for controlled 
unitary evolution generally includes frequencies close to the transition frequencies of the 
controlled energy levels (e.g., levels |1) and |2) for a two- level system). The radiation used 
for measurements includes components with frequencies close to the transition frequencies 
between the controlled levels and the ancilla levels, which moreover should be subject to 
decoherence upon decay to the lower energy levels. In general, measurements also can be 
performed through collisions between the system and electrons or atoms when the scattering 
cross-section depends on the initial state of the system. In this case the scattering data 
will provide information on the initial state of the system, thus realizing a measurement 
procedure. 

Measurements on a two-level system in an arbitrary orthonormal basis and can 
be realized using an ancilla system and inducing an interaction Hamiltonian between them, 
which generates a unitary evolution operator U$ such that for any vector \x) = <y\ip) + (3\ip') 
of the initial system one has L^|x)|l') = a|^)|l') + j3\ip') |2'), where and |2') are the 
energy levels of the ancilla system. The unitary operator can be chosen as U$ = I — 2P, 
where P = is the projector onto the one- dimensional subspace of the composite system 

spanned by the vector \(f>) = — \2'))/y/2. Then, nonselective measurement of the 

energy level populations of the ancilla system [i.e., measurement in the basis of |1') and |2')] 
realizes an indirect measurement of the initial system in the basis and and changes 
its state into p = \a\ 2 \i/))(i()\ + \(3\ 2 \^')(^'\- 

An indirect arbitrary von Neumann measurement on some quantum system can be exper- 
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imentally realized if the system can be coupled with another appropriate ancilla system (the 
ancilla can be identical to the initial system) and any unitary operator between these two 
systems can be implemented. For the case that the measured system is a two-level trapped 
ion, the ancilla could be another two-level trapped ion. Arbitrary unitary operators between 
the two trapped ions can be implemented using a sequence of at most three controlled-NOT 
(CNOT) gates and fifteen elementary one qubit gates 22J (i.e., single ion unitary oper- 
ations). Therefore experimental realizations of CNOT two-qubit gates [3] together with 
ability to realize arbitrary one-qubit unitary evolutions allows for generating any two-qubit 
unitary operator, in particular, the operator U$ from the preceding paragraph. Then the 
ability to measure the ancilla ion in the energy level basis makes arbitrary measurements on 
the initial ion practically possible. The detailed specification of such a scheme for practical 
laboratory realizations of optimal measurements from Sec. IIHI requires a separate study. 



III. CONTROL BY MEASUREMENT IN A TWO-LEVEL SYSTEM 

This section presents the analytical solution for maximizing the expectation value of any 
given target observable of a two-level system by optimized measurements. First, the case 
with neglect of system free evolution between the measurements is considered. After that the 
modification induced by including the free system dynamics is described. Any non-trivial 
observable Q of a two- level system is an operator in C 2 with the form Q = q\P\ + qiPi, 
where q\ and q 2 are eigenvalues and Pi and P 2 are the corresponding projectors such that 
P\Pi = and Pi + Pi = I. The observable Q is measurement-equivalent to the projector Pi 
(or to P2 = I — Pi)- Thus, any non-trivial observable of a two- level system is measurement- 
equivalent to a suitable projector and the problem of optimizing over the most general 
measured observables is equivalent to optimizing over measurements of only the projectors. 

Any density matrix of a two-level system can be represented as 

p = _[l + w .o-] 

where a = (01, 02, 03) = (<r x , a y , a z ) is the vector of Pauli matrices, w = (w\, 102,103) = 
(w x , w y , w z ) e M 3 is the Stokes vector, |w| < 1, and w • a = w x a x + w y a y + w z a z . Thus, the 
set of all density matrices of a two-level system can be identified with the unit ball in M 3 
(i.e., the Bloch sphere). Given a density matrix p, the components of its Stokes vector can 
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be calculated as Wi = Tr (p<Xi) for % — 1, 2, 3. 

Pure states correspond to projectors which can be represented as density matrices with 
Stokes vectors of unit norm, |w| = 1. Measuring a projector P transforms the initial density 
matrix p into the new density matrix pi = Aip(po) defined as 

Pl = P Po P + (I - P)p (I -P) = Po- [P, [P, Po]] 

Let a and a be the Stokes vectors characterizing the initial density matrix p and projector 
P, respectively (so that |a| = 1). Using the commutation of Pauli matrices [<7fc,crj] = 
2ieki m <7k, where e^ m is the Levi-Civita symbol, one gets 

1 1 3 

[P,Po] = -[(I + a-<7),(I + ao-<7)] = - ak(ao)iWk,0i] 

k,i=i 

i 3 ' 3 

= ^ e kimak(a )ia m = - ^(a x ao) m (J m 

k,l,m=l m=l 

where (a x a ) denotes the vector product of a and a . This gives 

px = i[I + (a + a x (a x a )) ■ a] 

Using the Lagrange formula a x (b x c) = b(a • c) — c(a • b) for the double vector product 
and noticing that |a| = 1 produces 

a + a x (a x a ) = a + a(a • a ) - a |a| 2 = a(a • a ) 

Therefore we finally have 

Pi = (a-a )(a-(j)], 

such that the Stokes vector of the density matrix after measuring P takes the form w x = 
a(a- a ). 

Consider a consecutive measurements of the projectors Pi, ... , P/v on the same physical 
system. Let a fe for k — 1, . . . , N be the Stokes vector characterizing projector P k . After the 
last measurement the density matrix will be 

p N := M Pn o M Pn _ 1 • • • o A^ Pl (p ) = ^t 1 + w n ■ <A 
Here the Stokes vector w^v has the form 

w N = a N (a N ■ a A r_i)(a A r_i • a/^) • • • (ai ■ a ) 

= a N cos at cos(pjv-i • • • c °s v?i|ao| (5) 
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where ipk is the angle between vectors afc_i and a&. 

The objective functional ([2]) may be rewritten as follows. The target Hermitian operator 
can be represented as = AoI+ A -ex, where Ao = Tr is a real number and A = Tr (Oct) G 

is a three-dimensional real vector. Using this representation produces Jjv[ a i> • • • ; a 7v] 
An ■ w y • A. The control parameters are the unit norm vectors a 1; . . . , a^r which determine 
Wat. Introducing the target vector Wt = A/|A| of unit norm, |wt| = 1, then the objective 
functional becomes Jn[&i, • • • , a v] = Ao + |A|(\Vjv • Wr). It is clear from this expression that 
maximizing the objective is equivalent to maximizing the scalar product w^r • w T . 

Let Aip be the angle between ao and the target vector wy and ifiN+i be the angle between 
vectors a^r and w T so that J^fJi 1 fi — ^V 9 - Here the inequality is used since vectors a fc may 
in general belong to different planes. The equality may hold if all vectors a , . . . , a^v and 
w T belong to the same plane. As a result, the objective functional may be expressed as 

Jjv[ a i> • • • > a iv] — Ao + |A||a | COS(y97v+l cos ip n . . . COS (fi 

The objective Jn is maximized i£ <f\ — (f2 — ■ ■ ■ — <fN — Vn+i = Atp/(N + 1) so that 
Y^fJi 1 fi = ^V 9 ! an d the maximal value of the objective is 

r Aw i^+i 

J^ = A + |A||a |[cos-^ I J (6) 

The corresponding optimal fc-th measured observables are those which are measurement- 
equivalent to the projector P£ pt = ~[I + a^ pt • a], where the vector a°, pt belongs to the plane 
formed by a and and is obtained by rotating the unit norm vector a /|a | by the angle 
kAip/ (N +1). Any such observable has the form 

QT = q k PT + Ul-PT) (7) 

where qt and % are real numbers with q^ ^ q\. It is not important for the control purposes 
here which observables from the equivalence class are chosen. In particular, the projector 
Pj° pt could be used as the measured observable. In general, if the free system evolution is 
neglected then all the vectors characterizing optimal observables must belong to the same 
plane formed by the Stokes vectors of the initial and final states. 

The Stokes vector of the system density matrix after measuring P£ pt , in the case that the 
total number of measurements is N, will be 



wjvfc = ao 



cos ■ 



N+l 
9 



(8) 



Thus, each optimal measurement rotates the Stokes vector of the density matrix by the 
angle A(p/(N + 1) and shortens its length by the factor cos[A<p/(iV + 1)]. 

The analysis above assumes that the free system evolution between the measurements 
could be neglected. This limit will be valid if the time between measurements At is suf- 
ficiently small such that AtEi <C 1 for the eigenvalues Ei, % = 1,2, . . . of H ; the limit is 
also valid if the system is close to degenerate, Ei ~ Ej, In order to go beyond this 

limiting case, now we will describe the modification induced by the free dynamics. Suppose 
that N measurements of observables Qi, ■ ■ ■ ,Qn are performed at the fixed time moments 
< t± < t2 < ■ ■ ■ < < T and between the measurements the system evolves with its 
time independent free Hamiltonian H . Then the density matrix p(T) at the target time 
T is given by the equation of the form (j3j) with the unitary evolution between (k — l)-th 
and k-th measurements U[t k _ lt t h ){p) = &^[—i(tk—tk-i)Ho]pexp[i(tk—tk-i)Ho]. The relation 
Mq(UpW) = UM Q (p)W with Q = WQU gives by induction p(T) = e~ iTH ° p N e iTH ° , where 
Pn = M.q n ° • • • °-M.Qi (po) is the density matrix evolved only under the measurements of the 
modified operators Qk = e ltkH °Qke~ ltkH ° . Therefore the objective function for a target op- 
erator in the case of including the free dynamics, Jn[Qi, ■ ■ ■ , Qn', 0] := Tr [p(T)0], equals 
the objective function Jn[Qi, ■ ■ ■ , Qn', ©] : = Tr [pat0] without the free dynamics with the 
modified target operator = e lTH °Qe~ lTH ° and measured observables Qk = e ltkH °Qke~ ltkH ° . 
The latter problem was completely solved above and the optimal measured observables 
are given by (EJ). If Q° k v are such optimal measured operators for the objective function 
Jn[Qi, ■ ■ ■ , Qn) 0] with neglected free dynamics, then the optimal measured operators for 
the objective function Jn[Qi, ■ ■ ■ , Qn) 0] are Q° k pt = e~ ltkHo Q k pt e ltkH " . This implies that, 
while all the vectors a^ p corresponding to the optimal operators Q° k pt belong to the same 
plane, this is not true for the vectors a£ pt corresponding to the operators Q^ pt . 

Thus, if the free evolution between the measurements is not important then, for arbitrary 
initial and target states, the vectors characterizing optimal measured observables must be- 
long to the same plane formed by the Stokes vectors of the initial and final density matrices. 
If the free evolution is relevant, then the optimal observables undergo additional unitary 
transformations with the generator H , which moves their corresponding vectors out of a 
plane. 



10 



IV. MEASUREMENT-INDUCED POPULATION TRANSFER 



Here we apply the general result to the problem of population transfer between or- 
thogonal states |1) and |2) of a two-level system. The initial density matrix is po = 
The target operator is the projector on the excited state, 6 = |2)(2|, which corresponds to 
Ao = 1/2 and A = l/2e 2 . In this case ao = — e z , wt = e z and therefore the angle between 
the initial and target vectors is A<p = ir (see Fig. [1]). The maximal population transfer to 
the excited level p<2,2{N) as a function of the number of performed measurements is given 
by ([6]) and has the form 



1 



7T 



N+l 



p 22 (A0 = -^l+^cos ¥TT j (9) 

The function p 22 (iV) for N < 50 is plotted in Fig. [2j 

Suppose that the goal is to transfer population p 22 = 1 — e to the excited level, where 
£< 1. The asymptotic number of optimal measurements necessary to meet this goal may 
be found as follows: 



p 22 (iV) = l-e=M 



7T 

cos ■ 



1 N+l 



2-2e 

N+l. 



N + l v ' N+l 
Notice that a small value of e requires N to be large. Therefore the Taylor expansion for 
cos(7r/(iV + 1)) may be used which gives the approximate relation 

7T 2 2e 



2(N+1) 2 N+l 

It then follows that the number of measurements necessary to transfer population p 22 = 1 — e 
to the excited level asymptotically behaves as iV w tt 2 /(4s), which is consistent with the 
behavior in Fig. [2] with N = 50 and e ~ 0.05. 

Figured] illustrates the population transfer by optimal measurements in a two- level system 
for N = 10. The left-hand plot corresponds to the case with a non-trivial free dynamics 
driven by the free Hamiltonian H = \<J Z . The 10 observables Qi, . . . ,Q±o are measured 
at the time moments t}. = Tk/(N + 1), where the target time is chosen as T = n, and 
characterized by the unit norm vectors 



-opt / irk \ / 7rk 



Tk \ ( Tk 

cos| ivTT>^ + sm (ivTl)^ 



n 



k= 1,...,N 
(10) 



shown on the left-hand plot. The smooth curve passing through the ends of these vectors 
represents the continuous family of the projectors characterizing the anti-Zeno effect in the 
limit of infinite number of optimal measurements, as described in the next section. 

The right-hand plot in Fig. [TJ illustrates the evolution of the system without free dynamics 
between the measurements. In this case ipi = ■ ■ ■ = ipu = and the optimal observable 
for the k-th measurement is characterized by the vector 

- cos f^M e - + sin f7F^) e - ( n ) 



a 



opt 
k 



N+1J * \N + 1 



which is obtained by rotating the initial vector a by the angle kipi. The system density 
matrix after the fc-th measurement has the Stokes vector = w 10j fc = (cos <£>i) fc a^ pt . The 
plot shows the length of decreasing after each measurement by the factor cos</?i. The 
relation between the two cases is that each vector a°, pt describing optimal /c-th measurement 
for the case with free dynamics is obtained by rotating a°, pt by the angle nk/10 around 
z-axis. Note that the ten vectors on the left plot of Fig. [TJ are the vectors a°, pt characterizing 
the ten optimal measured observables which are measurement-equivalent to the projectors, 
i.e., pure states, P£ pt = ~[I + a^ p • a] and therefore these vectors have unit norm. The ten 
vectors on the right plot characterize the system density matrix p = |[I+bfc-cr] after each 
measurement, which is a mixed state due to decoherence induced by the measurements, and 
therefore these vectors have norm less than one. The shortening of the vectors characterizing 
the system density matrix after each measurement for the case illustrated on the left plot 
will be the same, as for the case shown on the right plot, i.e., this shortening is by the factor 
cos(/?i after each measurement. 



Numerical simulations for some models in 



also suggest that the expression gives 



the estimate for the maximal population transfer between any pair of orthogonal levels in any 
multi-level quantum system. The complete analytical investigation of this problem remains 
open for a future research. 

V. RELATION WITH THE ANTI-ZENO EFFECT 

The maximal population transfer to the excited level by a finite number of measurements 
satisfies P22(N) < 1. Since lim [cos(7r/ri)] n = 1, one has lim p22^N) = 1, i.e., complete 

n^roo N—*oo 

population transfer is attained only in the limit of an infinite number of measurements 
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which describes the anti-Zeno effect. To show this behavior consider the projector E = 
|1)(1| = | [I + ao • a]. The goal is to steer the system from the initial state |1) at time 
t = into the target state |2) at time t = T. Define for each t e [0, T] the unitary operator 
U(t) = exp(ia y 7rt/2T). Then E(t) := U(t)EW (t) is the projector characterized by the 
Stokes vector w f = — cos(%t/T)e z + sin(7rt/T)e x , i.e., E{t) = + w t • a]. According to the 
anti-Zeno effect, continuous measurement of the projector E(t) steers the system at time 
t into the state characterized by vector w t . One has wt = e z , i.e., at the target time the 
system will be transferred into the state |2). 

The anti-Zeno effect is obtained in the limit of an infinite number of measurements when 
the interval between any two consecutive measurements tends to zero. Consider performing 
N measurements of the optimal observables Q° pt , . . . , , where Q°, pt = ~[I + a°, pt • a] and 
a°, p is defined by Eq. (11 II) . If the dynamics between the measurements is neglected then the 
state of the system after the k-th measurement will be characterized by the vector 

/ iik \ / nk \ ~ 

TT — r + sin — — - )e x 
\N + lJ VA^ + 1/ . 

Taking the limit as N, k — > oo such that the ratio kj (N + 1) = t/T is kept fixed produces 

lim wjv k = — cos(iit/T)e z + sm(irt/T)e x = w t 

Thus, the anti-Zeno effect is recovered in this limit. The corresponding evolution is char- 
acterized by the projector E(t) describing the rotation of the Stokes vector of the system 
density matrix in the same plane. In general, other evolutions exist which steer the system 
from the ground to the excited state with projectors E(t) which induce rotations out of a 
plane. Such projectors can be limits of optimal control by a finite number of measurements 
if the free evolution between the measurements is non-trivial. As an example of this sit- 
uation, the smooth curve on the left plot of Fig. [1] describes the optimal anti-Zeno effect 
characterized by the projector E(t) = ~ [I + w f • a] , where 

w t = lim a°, pt = — cos(nt/T)e z + sin(7rt/T)[cos(7rt/T)e :r + sm(irt/T)e y ] 

Here the vectors a°, pt defined by (flQj) characterize the optimal measurements for the example 
with the free evolution considered in Sec. HVl the target time is chosen as T = tt, and the 
limit is taken with fixed ratio k/ (N + 1) = t/T. 
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VI. CONCLUSIONS 



In this paper control of quantum systems by non-selective measurements is considered. 
The capabilities of optimized measurements for control of a two-level system are explicitly 
investigated. The optimal observables and maximal expectation value of any target observ- 
able are analytically found given any initial system density matrix and any fixed number 
of performed measurements, thus providing a complete analytical solution for control by 
measurements in a two-level system. For any given number of measurements the degree 
of control, i.e., the maximum value of the objective, is found. The relation between the 
optimal measurements and the quantum anti-Zeno effect is established. Looking ahead, the 
ultimate goal will be specification of laboratory protocols to make the procedure of control 
by measurement practical for realistic systems. 
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X 

FIG. 1: The plots illustrate evolution of a two- level system during the optimal measurement- 
induced population transfer from the ground to the excited state. The vectors ao = —e z and 
wt = e z are the Stokes vectors of the initial and target states. The total number of measurements is 
N = 10. The left plot shows the ten unit norm vectors describing the optimal measured observables 
for a case with non-trivial free system dynamics. The smooth curve passing through the ends of 
these vectors describes the projectors for the anti-Zeno effect in the limit of an infinite number 
of optimal measurements. The right plot shows the Stokes vector = w^jt, defined by O, of 
the system state after k optimal measurements with neglected free dynamics. Each measurement 
rotates the preceding vector by the angle (p% = anticlockwise and shortens its length by the 
factor cos</3i. 
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1 




Number of measurements, N 

FIG. 2: Maximal measurement-induced population transfer from the ground to the excited state of 
a two-level system as a function of the number of optimal measurements. The maximal population 
P22(N) approaches 1 as N — > oo. 
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